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ABSTRACT 


In view of the recent interest in mathematical models 
simulating food chains, the present dissertation is devoted to a 
study of some aspects of predator-prey food chains by means of 
deterministic mathematical models. These aspects include: 
boundedness of species numbers, stability of equilibria (both 
local and global), persistence and survivability, alternating 
equilibrium numbers phenomenon, perturbations (autonomous and 
periodic), effect of selection on the abundance of species, and, 
relationship between stability and complexity. My main results 
are Theorem 4.4.1 and Theorem 4.6.1. There I show that for Lotka- 
Volterra food chains of any length n with a carrying capacity 
in the first species, the positive equilibrium, if it exists, 
will automatically be globally stable. Hence there cannot be 
any non-trivial limit sets such as periodic orbits or recurrent 
motion. Moreover, if the parameters are varied in such a way as 
to make the positive equilibrium disappear, then at the instant 
of disappearance of this positive equilibrium it coincides with 
the equilibrium of one dimension lower, making the latter globally 


stable and so on. 
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CHAPTER I 


INTRODUCTION 


§1.1 Background and Survey of Literature 


Recently there has been considerable interest in mathematical 
models simulating predator-prey food chains. In the case of food chains 
of arbitrary length, the most general Kolmogorov-type model is looked 
at by Rescigno and Jones [1] who give some hypothesis and geometrical 
interpretation of the equilibrium leading to oscillatory phenomena, 
Saunders and Bazin [2] consider the stability and the effect of enrich- 
ment on the stability of a food chain in a chemostat using an inter- 
mediate-type model. More recently, Gard and Hallam [3] give sufficient 
conditions for the persistence and non-persistence of Lotka-Volterra 
type food chains. For food chains of length three, Rosenzweig [4] and 
Wollkind [5], using graphical methods, analyze the neighborhood stability 
of the possible equilibria of a general food chain model simulating the 
plant-herbivore-carnivore dynamics. Freedman and Waltman [6] examine 
the survivability of all species of an intermediate type model. MHausrath 
[7] considers a perturbed Lotka-Volterra type model in which the highest 
trophic level is not soley limited by the middle level and he shows the 
existence of a perturbed asymptotically stable equilibrium. More 
recently, Lin and Kahn [8] discuss the existence of a bifurcation of 
equilibria and calculate to the lowest approximation the periodic 
solution around one of the unstable equilibria of food chains having 


Lotka-Volterra type and Holling type of predation. Food chain models 
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where both biomass and energy have been taken into consideration can be 
found in May [9] and Hirata and Fukao [10]. 

All the food chain models considered by the above authors can 
be included in the following ecocommunity model of n variables (c.f., 


Levins [11]): 


x' = £ (ky Xp oe eo : Cj 2C,5+-+) 


Ke- 


xi = £ A(X] Xo 0002 9X : Cj 5C,5++-) 


No = 


Gi.41.1) 
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d ; 
where ' — and t denotes time. 


~ dt 
That is, the growth rate of Xs is some function (called the growth 
function) f. of all the variables Xp oXqo+00 9X) in the system and a 
set of parameters Cy 2Cos-es 3 

Here, for simplicity, we shall assume the variables x, 
(i = 1,2,...,n) is the abundance (number or biomass) of the eel. 
species. Generally speaking, the variables can also represent the energy 
content of the species, the level of resources used by the species or 
some factor produced by a species such as a toxic by-product or a pre- 


dation pressure. The parameters C,,C are of several kinds. Some 


C22°°° 
are environmental parameters beyond the influence of the species such as 

temperature, rainfall or vegetation density. Some depend on the genotype 
of a species, on such inheritable traits as its fecudity or heat tolerence. 


Still others are dependent on more than one species. For instance, the 


rate of predation of es on xX, depends on the biology of both. 
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By using (1.1.1) as our starting point, we have made the 
implicit assumption that the size of each population can be measured by 
real numbers, that no random fluctuations are present and that the age 
structure composition of the population is constant or irrelevant. (For 
a justification of using deterministic models such as (1.1.1), instead 
of the probably more realistic stochastic models, one may refer to 
Smecie ty h2 |). 

Let C1 2Cy5--- e P (the parameter space of interest) be fixed. 


Then (1.1.1) becomes 
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We shall assume that 


p n ne Thee ; 
(i) foofjs---ot, : Ry — R, where Ry =O ee Rares x, 00m WV eis cents 


are continuous and sufficiently smooth to guarantee that the initial 
value problem for (1.1.2) with initial condition €«€ Re will have 

a unique solution. Note that we only require the solution to be 
defined on some interval [0,T) where T>OQO may be <@, 


(Standard sufficiency theorems for existence and uniqueness are 


given in Coppel [13]). 


(ii) The hyperplanes H. = {xe Ry ce li 0} (Anes tJ ee. 


are invariant sets of (1.1.2). 
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This) 1s7itrue., sorsexample,; if (1.1.2) is of the form 


x = x) B(x) X50++ 29k) 
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Clearly, this assumption is biologically reasonable. 

Since the state space of ecological interest is the invariant 
set Be instead of ne notions such as equilibrium points, stability, 
asymptotic stability, global stability, periodic solution, recurrent 
Motion, etc. refer to solutions of (1.1.2) in Ee . Their mathematical 
definitions are given in Bhatia and Szego [14]. Discussion of some of 
these concepts within the ecological setting can be found in May [15] 
and Lewontin [16]. 


Following Gard and Hallam [3], (1.1.2) is said to be persistent 
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where [0,T) is the maximal interval of existence of x(t). In 
a. like manner, (1.1.2) is said to survive if each solution x(t) of 
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where [0, T) is as before. 


Remark: An alternate, but related, definition of persistence is 


given in McGehee and Armstrong [17]. 


Our main purpose in this dissertation is to study the stability, 
persistence and the effect of evolution on predator-prey food chains 
with special emphasis on Lotka-Volterra food chains and food chains of 


length four. 
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§1.2 Mathematical Preliminaries 


In this section we shall collect some results from ordinary 
differential equations and matrix theory which we shall have occasion 
to use in the following chapters. Definitions of the undefined terms 


and further details can be found in the references given. 


(1) Dulac's Criterion: This criterion will be used in §4.6 to show 
the non-existence of a periodic orbit for Lotka-Volterra food chains 


of length two with a carrying capacity in the first species. 


Statement of Criterion: Let 


ee FG Px) 
(1.2.1) neha ae: 


1 ~ £9 (&,> %») 


be a ct dynamic system and G a simply connected subregion of the 
domain of derinition of system. (1.2.1). Let B(x), xX») bea ct 


function defined in G such that the function 


5 5 
ax, (B £,) + ax, (B f,) 


does not change sign. Then G contains no non-trivial periodic orbit 


rays YORE YI 


Reference: Andronov, Leontovich, Gordon and Maier [18], p. 205. 


(2) Routh-Hurwicz Criterion: This criterion will be used in §3.2 and 
Chapter VI to obtain conditions for the asymptotic stability of 
the positive equilibrium. 

Statement of Criterion: The table below lists necessary and sufficient 

conditions for every root of the real polynomial in the left column to 


have negative real part. 
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Reference: Coppel [13], p. 158. 


(3) Schwarz's Theorem: We shall use this theorem in §2.2 to study the 
asymptotic stability of the positive equilibrium of food chains 


whose variational matrix has a particular form. 


Statement of Theorem: If a n Xn matrix A= fa, 5] can be transformed 


using a similarity transformation to the form 


(so called Schwarz matrix), where the unspecified entries of B are all 
zero, then the number of eigenvalues of A which have positive real 


parts is equal to the number of positive members in the sequence 
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does not change its eigenvalues. 


(ii) Eigenvalues of -A are the negative of those of A. 
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Reference: Schwarz [19] and Barnett and Storey [41], p. 52. 
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(4) Qualitative Stability Criterion: This criterion will be used in 


§2.2 and 83.2 to obtain necessary and sufficient conditions for 


the qualitative stability of the variational matrix. 


Statement of Criterion: A nxn matrix A= la,,! is qualitatively 


stable (i.e., matrices which are sign equivalent to A have all their 


eigenvalues with negative real parts) iff 
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Reference: May [15], p. 71 and Quirk and Ruppert [20]. 


(5) Barbashin-Krasovskii's Theorem: We shall use this theorem in 
§4.6 to analyses the global stability of all the possible equilibra 
of Lotka-Volterra food chains of any length with a carrying capacity 


in the first species. 
Statement of Theorem: Consider the vector differential equation 


ree se) 


with F: R = Re continuous and F(0) = 0. Let there exist an 


infinitely large positive definite function V_ such that the trajectory 


derivative V' of V_ satisfies 
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vi <0 outside M 


v' <0 on M 
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where M (cR ) contains no entire trajectories except the equilibrium 


point x=0Q0. Then x = 0 is stable in the large. 


Remark. This theorem can be trivially adapted to the case when the 


: F n 
State space is not necessarily R. 


Reference: Barbashin and Krasovskii [21] and Barbashin [22], p. 49. 
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§1.3 Outline of Thesis 


In Chapter II, we discuss general food chains, in particular, 
we give some ecologically reasonable restrictions on the f,'s in 
(1.1.2) so that (1.1.2) can be used as an effective model for predator- 
prey food chains. Some general theorems on the stability and surviv- 
ability of general food chains are also given. Chapter III is concerned 
with intermediate food chains. Boundedness of solution is proved in 
§3.1 whereas conditions for existence and stability of the equilibrium 
points of food chains of length four is stated in §3.2. In Chapter IV, 
we study the global stability of Lotka-Volterra food chains using a 
Liapunov function. We also illustrate the close relationship between 
stability and survivability in this particular case. Some easily 
obtained perturbation results are also given. Chapter V deals with the 
effect of selection on the equilibrium numbers of a food chain. Besides 
working out the details for food chains of length n = 3, 4, we also 
make some observations for food chains of an arbitrary length n. The 
three examples given in Chapter VI are used to illustrate the model 
possibilities in the controversal topic ‘stability verses complexity’. 
Finally in Chapter VII, we make a brief summary of the results that we 


have obtained and discuss their significance. 
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CHAPTER II 


GENERAL FOOD CHAINS 


§2.1 Assumptions 
Following Rescigno and Jones [1], we shall give a set of 
assumptions so that (1.1.2) can be used effectively to describe the 
dynamics of a 'straight' predator-prey food chain. 
Verbally, the assumptions are: 

(i) If predator and prey are both small in number, the prey will 
multiply (being unhindered by the small number of predators). 

(ii) The prey cannot multiply over a certain population size (its 
carrying capacity), even in the absence of predators, due to limitation 
of its own food supply. 

(iii) If the predator is too abundant, the prey cannot multiply. 

(iv) Multiplication of the prey is slowed by an increase in the 
number of predators. 

(v) For a constant predator/prey ratio, the multiplication of the 
prey (resp. the predator) is slowed (resp. accelerated) by an increase 
of the population size because the predator-prey encounters are more 
frequent. 

(vi) If there are not enough prey, the predator cannot multiply. 

(vii) Multiplication of the predator is slowed by an increase in 


its own numbers. 


We now translate assumptions (i)-(vii) into mathematical restrictions 
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(1) Since we are only concerened with straight chains, i.e., xX, Preys 


on xX 


1 *3 Preys On Xy5---, and x, Preys on x 1, we have, 
£5 (x) sXyo+++5%) = £ (x) >) 
£ (x) Xp 00+ +9%) = £ A (&) 2X5 9%) 
£ (x) Xy5+++9%)) = £ (Xo »X39X,) 
aed) 
£ ni (*%12% 22° + * 3%) : NES ie os oe 
f (8p oo 0+e 2%) = tO Sy-12* 
of of of 
CZ.) — <0 == coal Oi 3 aa <0 (€ because of (iv) ). 
2 3 n 
of of of of 
(3) = —J0', = MO SHE. el SU oa) 
2 |x,=0 ae <a n-l |x =0 S 
3 4 n 
(because of (vii)). 
of, of 
Note that, for example, we require i —. <0) ©instead sor <a) 
Xo x4=0 OX, 


because for X, > 0, increase in prey may satiate predators and thereby 


benefit the prey. 


(4) £ (0,0) > 0 (because of (i)) 


Hover Oy csuch jthat £(A;0) 0) (because of (ii)). 
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(6) mR 0, mR SC vcr acn ys oR <1i05 mR <0 where Ri is the 
i Dd n-1 n 
al 
direction of the line a constant in the X1-X5 plane; 
2 
ws 
R, is the direction of the line or ie constant, x, = constant 
S) 
in the Xj—Xy-Xy space; 
R-1 is the direction of the line = constant; 
n 
Sais constant in the X27 ¥n- 1 *n space; 
*-1 
and RY is the direction of the line = = constant in the 
n 
a a plane 
(because of (v)). 
Co ik Si> Sy» rats aa) > 0: such, that 
f,(S,;5 0; 0) =0 
£(S,, O50) = .0 
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ky : [s,, oo) > R 
ky : [S,s co) > RQ 
Le eee a 
such that 
£,(x,> Ox k3(x,)) = 0 V x aie 
£.(x55 On k) Gx) = 0 Vv X» 2S, 
£1 p-1? o; k(x)? 40 y ae? 420 


(because of (iii) and (i)). 


Remark: Although the assumptionsgiven above are ecologically reason- 
able, they do not include all the food chain models studied in the 
literature; for example, the Latka-Voterra food chain (see 84.1) is 


one important exception. 


y i Jb ig 
hatte 
Let met) 
i] 
j 
ay to « 1 : 4 y 
- . A t 
i" ie 3 Bigd 
? a : 
| 
J F “% ’ 
a {7 
t 
‘ 
Z \ 9 ace7 - ’ a i — | 
a i # ¥ . } | ’ b -«* ! 
fas" “a2 | Bs >" Weg ere 
} Fe 
4‘ y} is + low 
—inibey wi ieo toi 65 Sal) eyes stay i ube) Lars ee: ais he ‘al 


+ 


mE TE AD S25) pabse bagi tanks emit ied orate Shad 


Bf Bo We5Suis. 6) 8a ALAA fod Ot Paes an aoe i i ti, f 
‘ ~ 4 ~~ 7 


d s 4 4 a} 
page ne naa 


»@ ad 
, oT ny 
pie = 
7 : ene j 
pei oo 
q quy : 7 a 
- « oP - ’ oun Ae 
Z : 7 r aA i 
= fl . . _s _ 
7 * rt — i 
r a , 
Gan 4 


Baie & 


Cee cak) 


A Stability Theorem 


If assumption (1) of §2.1 holds, then (1.1.2) becomes 


' = 
x) = £, Gx) »x>) 
' = 
X5 £4 (x1 Xo 2X3) 
' = 
xy = £,(x).x3.%,) 
' = 
Seat Ent p22? *p-1 *n? 
’ = 
an £ Oy >%)) 


Suppose (2.2.1) has a positive equilibrium x*. 
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The following theorem gives necessary and sufficient conditions 


for M to be qualitatively stable. 


Theorem 2.2.1: M is qualitatively stable iff 


Mee he ae 
fife Wy he Soe Pee a such that m,, < 0 
and 
det M # 0 
Proof: It follows immediately from the qualitative stability criterion 


Stated in 31,2. (4). O 


Remark: Using induction (c.f., 85.3), one can show that for 


Mi? M5» » mi sO) then 


8) if ny eke leven 


Ciy) det M 
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Corollary 2.2.2: If m4 <0 and Mot Maa 7 eee SS 0 then M 
is qualitatively stable, and hence x* is asymptotically stable. 
Proof: It will be proved in 85.3 that under the above assumptions 
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The corollary then follows from the above theorem. JU 
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§2.3 General Theorems on Persistence and Survivability 


Of the numerous ecosystem considerations, probably none affects 
ecosytem structure as drastically as the extinction of a biological 
species (e.g., Paine [23]). As demonstrated by public concern (Holden 
[24]) extinction of a biotic component of a food chain is topical as 
well as ecologically fundamental. It is the purpose of this section to 
State and prove some sufficiency theorems on persistence and surviv- 


ability which are applicable to general food webs. 


Derini tion: =|Let.-"w Ry, + R be continuous on RQ,» e eyy kaye (Ri), 


and w(0) QO. Suppose also that the interval of definition of every 


solution y(t) of the differential equation 


(esol) y' =wy) 


with y(0O) « R, can be extended to the whole semi axis [0,~). Then 


(i) (2.3.1) is said to be of persistent type if every solution 
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Example: If w(y) = ay, then (2.3.1) is of persistence type whenever 


a > 0, and is of extinction type whenever a < Q. 
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y Mis w Gre, y(O) « Ba, 


has a unique solution, lim y(t) always exists (may be +), 
[ep aeae 


(ii) There are, of course, functions w for which (2.3.1) is 


neither of persistent nor of extinction type. 


eee nN . 
Definition: Let 0: Ry > R, be continuous on Ry and ct on 


n 
n - a ele) 
Ine (R,) and let o0 (x, 5+++5%,) = i pore f(x. 


-,x.) be the 
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PrvasectOory derivative.of p along (1.1.2). Then 6 is said to be a 


(i) persistence function for (1.1.2) if 


P(X, 5+-+5%)) >.0 whenever xe 0 
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where w is of peristence type. 
Giipy “extinction function for (1.1.2) af 
P(x, ,+++5%,) = 0) only if x, an), 


for some i= 1,...,n3; and, o' satisfies the inequality p' < w(p) 


where w is of extinction type. 
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inequality p' >w(p) where w is of persistence type. 


Theorem 2.3.1: (i) If there exists a persistence function for (1.1.2) 
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tay nou-persiatent (i.e.,  Gl.1.2) is not persistent); 
(iti )eitetberesexictspaucurvival tinction, fors(1, 1.2) ,9then, (1.1.2) 


survives. 


Remark: Parts (i) and (ii) of the above theorem are due to Gard and 


Hatlam [3]. 


Proof of Theorem: Since the proofsof all three parts (i), (ii) and (iii) 
are very similar, we shall only prove (iii). 
Suppose (1.1.2) does not survive. Then there exists a solution 


Mit) erok #GL. a2) with» -x(O) € Int (RL) such that 


fede leec Loca) WEE elim x, (t) = 0 
eT 


where [0,T) is the maximal interval of definition of x(t), i.e., 


there is a sequence {t,} with t, >T , such that 
x, (t,) > 0 as k>o . 


Since op is of persistent type, p (x(t, )) > Oi asamices eo" 50) OTethe 
other hand, the differential inequality 

py == we) 
implies o(x(t)) > y(t) where y(t) is the solution of the initial 


value problem 
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This contradicts the fact that w is of persistent type. J 
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in the case of Liapunov functions, it is in general not easy to 


construct a suitable persistent or survival function for (1.1.2). 
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CHAPTER III 


INTERMEDIATE FOOD CHAINS 


§3.1  Boundedness of Solutions 


In this chapter, we shall discuss intermediate type food chain 
models studied by Freedman and Waltman [6] (see also Freedman [25]). 


They are of the form 
x7 = *18(x)) - x54, (x) 
*5 = Xy(-dy + C54, (x1)) - X345 (x5) 
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Here the d. Ss \(i-#?2eo82 0) Sete athetdeath rate, of the i species if 
left alone, g is the self-generating function for the primary producer 
X10 Py (i = 1,...,n-1) the predation function of the ae species by 
the Ca)" species and C, (i = 2,...,n) is the efficency of pred- 
ation. 
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§3.2 Existence and Stability of Equilibria for n= 4 


In this section we shall find conditions for the existence of 
equilibria of the intermediate food chain model (3.1.1) when n = 4 
(which is the plant-herbivore-carnivore-secondary carnivore situation). 
Also, we shall discuss the asymptotic stability and the instability of 
these equilibria. 
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To study the stability of these equilibria, we consider the 
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Hence, Ey, Willebeuwasyuptoticallvyestable 46 (3.2.3), (o22.4))., 945: 259) 


and (3.2.6) hold. We also have the sharper result 
Theorem 3.2.1: M is qualitatively stable iff 


and M3 <0 
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Proot: Note that M)4™59 a Mj 9M51 is equivalent to the condition that 
0 not be a characteristic root of M (i.e., det M # 0), which is 
puaranteed by the above conditions (since, in that case M,4™)5 EN 


The remainder of the proof follows from Theorem 2.2.1. O 
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CHAPTER IV 
LOTKA-VOLTERRA FOOD CHAINS WITH OR 


WITHOUT CARRYING CAPACITY 


84.1 Introduction 


Lotka-Volterra dynamics has been extensively studied by many 
authors. Besides the classic works of Volterra [26] and Lotka [27], 
one may also refer to D'Ancona [28], Kerner [29], Goel, Maitra and 
Montroll [30], Balinese and Richardson {3.|, and Goh [32]. in this 


chapter, we shall consider Lotka-Volterra food chains of the form 


' — _ _ 
X, = X3(-agq + a39%) ~ 239%3 — 454%,) 
Ge a) 
! ae a = a 
steal is x16 fn-1 On net no2eneo oat n=1onel n= n*p? 
x x, (-aig Bs a Ante i an n x,) 
where all the specified a8 are > QO, except ai anor see an 


which are > 0. 

Our purpose is to study the stability of the equilibria of 
(4.1.1) and to improve some of the persistence results given in Gard 
and Hallam [3] to survivability. 


It is clear that the two assumptions given in $1.1 hold 
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Por (4. App orine 


n : : ; 
V Xo € Ry » there exists a unique solution x(t) of 


(4.1.1) such that x(0) = Xo 
and that H, (iverd ci.) tare, anvartante setsrob, Cael 


Moreover, we have 


Theorem 4.1.1: Jive a1? OF * thenvevery solution of -(4.d.1)) with 


om ee oh NM. 
initial condition «¢ R, is bounded. 


Proot : (Similar to Theorem 3.1.1). U 
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§4.2 A Liapunov Function 


Lemma 4.2.1: System (4.1.1) can be rewritten as 


n 
(42.1) x' = x,(e, - A.x, + = ) cGy BS) (ial 
al 2b ak B. j=1 1j J 
where €) = 419° €2 = 74299 see ta 5 
ds = a., (4 2 “jane Be >) (i = 1, 


and Fora) is anti-symmetric. 


PrOoLs From the first equation of (4.1.1), choose 


1 = = _ 
B = 1 and a aio 


From the second equation of (4.1.1), choose 
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we have the desired result. OQ 


Suppose there exists a positive equilibrium 


x* = (x#; xf Sets x*) foreG4. 2:8). hence alsovftor (401.1). - Define 


n x 


a — Ct oe * —_—_ 
V(xysXpo0e00%,) = 2 By Oxy — xR - xf In xe? 
i=l 
(for XyoXooe009%, ? 0). 
Lemma 4.2.2: (i) V(x*) = 0 and 


VG), 4240 for x # x* 


(ii). The ftrajectory derivative Vv’ tof. V salong (4.52.1) dis 


given by 
e 2 
! ees = xk 
V (KX) oXy0+++5X,) A ds Bs Gx, x*) 
i=1 
(iii) V(x) >* as x tends to any of the invariant hyper- 
planes H. (feo eo ay OF ase hex ieee 
Proof: Consider the function f : (0,~) > R defined by 
x 
f(x) =x-a-atln F (x 20) 


where a> 0..is fixed. 
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Therefore, f attains its minimum at x = a; and (i) then 


follows. For (ii), 
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City is clear, «.O 


Theorem 4.2.3: If a positive equilibrium x* exists for (4.1.1), then 
Gi), =x*— fs stable, 


(ii) solutions for system (4.1.1) with initial condition 


é int (Ri) are bounded, 
(iii) system (4.1.1) survives. 


PrOoL: (ii) and (iii) follow from Lemma 4.2.2. (i) follows from 


hema 4.222 and Hale [33], p.. 293, Theorem 1.1. 


Theorem 4.2.4: Suppose there exists a positive equilibrium x* for 
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(i) x* is globally asymptotically stable (over Int (R,)) 


(44) ihe varshational matrix M of (4.151) at x* is qualitatively 


stable. 


Proof: (i) follows immediately from Lemma 4.2.2 and Hale [33], p. 297, 


Gorolmary 1.2... ror (11), observe: that 
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Therefore by Theorem 2.2.1 and the remark following it, M is 


qualitatively stable. JQ 


Remark: Using determinants, one can write down conditions for (4.1.1) 


to have a positive equilibrium. 
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§4.3 Persistence and Survivability 
In this section, we shall consider the particular case of 


(4.1.1) when Aono = oes Ba Gaye ric ear: 


11 *1 7 412 *»? 


Hp = X_ (agg t+ any ¥y — An3 3 
tice fr ma 
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(4.3.1) 
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Theorem 4.3.1: 
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Then there exists a positive equilibrium x* for (4.3.1) iff Ug > 0, 
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a =a and = 0 if n= 2m 


01 11 Fomt1 0 


(iii) The positive equilibrium x* in(i) and (ii), if it exists, 


is unique. 


Proof: The proof of (ii), which is similar to that of (1), will be 
given in §4.4. (iii) follows from the proofs of (i) and (ii). Consider 
G). x* = (x#, ne eT Ae x*) is a positive equilibrium for (4.3.1) iff 
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aie ee Mee Take are >Q. On the other hand 
if 3 n=-1 
a a 
ae ,* n-1 Q i n-l n og 
aor Po-1 n=2 fos lsneoeee 
a a 
xk = _ 30 wes x 
32 #32 
all X59 « »x* 9X8 ares 0) ff xF > 0 Thus x* exists 
Za0u. lo-calculate x* we note that 
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Corotiarn “4,302: Let n, be even in (4:3.1) with 4447 Oe e 
oO then 


ae 
(i) a positive equilibrium x* of (4.3.1) exists, is unique 
and is stable, 
Cay tevery solution of (4.5.1) with initial condition “e Int (Ry) 
is bounded, 


(iii) system (4.3.1) survives, 
(iv) each of the hypersurfaces 
Nes = {x e Int (Ry) SV) eed where c>0O 
contains a recurrent motion for (4.3.1). 


Proof: (i), (ii) and (iii) follow from Theorem 4.2.3 and Theorem 4.3.1. 


i = =... = 0, Lemma 4.2.2, (V' =0, so thar V 
Since ary ano a . by “d 
is a non-empty compact invariant set for (4.3.1). (iv) then follows 


from a classical theorem of Birkhoff (see Nemytskii and Stepanov [34], 


p. 374, Theorem 7.02 and p. 375, Theorem 7.06.) U 


Remark: The variational matrix of (4.3.1) at x* with 4417 0 is 
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of the form 


0 {2 
95) é api: 
30 


By the same transformation as in the alternate proof of Corollary 2.2.2, 


M can be transformated into the Schwarz matrix 
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Therefore all the eigenvalues of M are pure imaginary. 
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bounded. 
Cit) a oait) a 220 .<8then 


(a) a positive equilibrium of (4.3.1) exists, is unique and 


is asymptotically stable, 


(ob) (4.3. 1) survives. 


EroOoL: We have only to show x* is asymptotically stable. The 
rest follows from Theorem 4.1.1, Theorem 4.2.3, and Theorem 4.3.1. 


Now the variational matrix of (4.3.1) at x* is given by 
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Mi) = yg ~ 28,1 x* - a) xt = —a, xf <0 
“apy casi es ieee 
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m = x* > 0 


Therefore Corollary 2.2.2 applies. U 


Remark: Gard and Hallam [3] proved that when ar4 or 0 pe (4eS AL) eis 
non-persistent if wu < 0. We shall show in 84.6 that (4.3.1) is 

non-persistent if yu < 0. They have also proved that when a = 0 
and n is even, then (4.3.1) is non-persistent if Uo <0. We now 


show that (4.3.1) is non-persistent even when Ug = QO. 


Theorem 4.3.4: Let n_ be even and 4447 OA ni a eee lee Uy = 0% 


then (4.3.1) is non-persistent. 
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ight 10 
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The variational matrix of (4.3.1) at x is given by 
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since x = 0. 
a a : , 
Now a0 Be x-1 is the eigenvalue of M in the 
a 
x direction and it can be made < 0 if we choose (0 <) x aS nO 
So nel 
*n0 
Therefore, for those points of L with x-1< ea each of them 
n n-l 


will have an asymptotically stable manifold intersecting Int CED). 


Hence (4.3.1) is non-persistent. oO 


Remark: In the case n is odd and 4447 0, Gard and Hallam [3] 


showed that Ug > 0 => (4.3.1) persistent and Ug <eQ ee sr(4,S 01) 


non-persistent where the definition of Uo is exactly the same as 
that given in Theorem 4.3.1, except that m= =. 
For Uy = 0, Freedman and Waltman [6] showed that when 


n = 3, every solution of (4.3.1) with intital condition e¢ Int (R?) 

is periodic, and hence (4.3.1) is persistent. The situation is unknown 
when n is arbitrary. However, it can be shown that when Ug = QO; 
there is) an infinite bine of equilibria of (4.3.1) in Int (RL) 


which is initiated from an equilibrium x of (4.3.1) in 
nh 
: wears 0 d == () 
{x € Ry: X1> Xo» es > an x } 


Random fluctuations can move the system from one equilibrium point to 
another. Eventually, these fluctuations will move the system to a 


neighborhood of x, and x, can be considered as extinct. 
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84.4 Equilibria for (1) 


In the following three sections, we shall be concerned with 


eruaetor (4.3, l)yawith ai 2.0, i:e%; 


a 


t = — 

9 = Ky agg + a57%) ~ 499%3) 

Be are SCR meee Shee ete | vee Qh ine 0) 

(1 ) oh aki pmnccaisl ane tay 
n 

x. =a (-a +a x 25) x ) 
n-1 n-1 n-l 0 n-l n-2° n-2 n-l nn 

x' = a 


n x, (ang n era la et 


where all the specified a are > 0. Consider the system of equations 


1 (49 ~ 441%] 7 242%2) = 0 


Ky (-aggq + a54%1 ~ 493%3) = 9 
C40 i 51) i 
x, 6-400 if a noha 0 
Since for Xp Xoo sees x aD 5 
= = = = => = 
x, = Q => X5 0 arr x OF; 


in (4.4.1), the only possible equilibria of (1) in R are of the 


form 


10 xoi a 
: Plaie Sa ged ang iat a 


LA ewer repye 
Ris" PD dob tel 


phe 
Nt { Aa. 
f a f 4 ae ; { a ey hae as 
oy L mt h , 
= ait 7 
: . 2) e“28 Me eee 
NS ss 


sp a ant Ob aaa 


A a 
o D 
| fells . é 
Ws Hh, ED) See) nradd tlope aie reeeaes 
/ ibs i 
A re 1 
& : - 
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Ey (Onm40' 50: Apt 18), 
By” a5 (F910, 0, me ae 0,0) 
ed a (*, de 0, - > 0; 0) 
(n) k ok ok x 
ey = ( cy 3 c) ’ 3 0) 
CUS ED i ae kok 
EY stl CAP EES tye 9] 


(n) 


where * denotes a positive number. Of course, EO is the positive 
equilibrium of (I). 


Our purpose in this section is to derive necessary and 


sufficient conditions for the existence of the equilibria 


(n) (n) (n) 
0°? EL Bia Paces EY of (I) 


E 
and to see what are their coordinates if they exist. 


th 


Notation: For xe a let Pp, (x) denote the i coordinate of 


x. (This was denoted by x). 


We separate into two cases. 


Case 1: n is even, n= 2m Say. 


Let Ez = x* = (x¥,+.+,x4) = (xf, BAe xt): Then since 
x* x are all > O and they satisfy 


_* ge 


. “oo an -. +e gel eai a es soakb yal 
he a? a4) 

aur Oa yal ‘eyise vr at endian wa 

ov b gap a ard a nah Sixes Sia sige eiyiethaos 


C33) ta} uh 
, or all q ve Wy , 


in oe 
em -ofnty EE nls git6 robs “i pers Pun oid 


ari? S5dreb! "CD fe peer ia ce 


a Pad re hie ie 
i i “i wa 


Ge ee, a sea 
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Rone) Oo fms) 2m? me ee noe 


-a =, Q 


sth 2 oe bees 


we can solve for them in the order 


* * * * * * * 
EAST OOS, eae Rm A I I Fe 
* * 
PO ame 2 tr 2 em 


Using a computation similar to the proof of Theorem 4.3.1 (i) we see 


that 


xT 5 Meee etavess x3 1 are all >0O 
* i i * 
x5 x7» Satie xB 2? x3 will ails a ba a x5 oO. 
Moreover, 
m a i-l a 
? 5 ay 
Cee. Pog ein a se &. it I oe (1 <k<m) 
oe f=k “21 21-1 j=k “23 25-1 
a a m aa 0 Oo DWE ea Spt 
1 240 2} 2jt1 
(4.4.3) p (EM) = xg = 10 “Woy _ 2i0 9, 23 25+ 


O10) Bl nte O21 el 24 5a 


and 


ede 


Fa x83 a a 

74 “hy q1 eh tn 

Ae Vg a he I ee oe 
ex + rag th ee. ink 


ey AG Gs ee 


ae 


ih 


Kits.. 
(4.4.4) peete = xh, =( lies cee ie see aN 
wi) Cee tera el] 

m Ler ae : 

i=l Had» 2h 2y 1 

tg 7 peer 

ON Ou ae eens ees) 

oa j=l 72541 25 

Convention: ) =0, TM =1 (6 denotes the empty set). 
> > 


Case 2: n is odd, n = 2mtl say. 


Similar to Case 1, we obtain 


Meare al ae Z 
(4.4.5) eo) = y a2tt0s II ee ial be (1 <i let .< m) 


P eas 
zt isk *2i+1 24 j=k+1 724-1 25-2 


a a Uy Br: Wee i a8 
ee ae My meen ee ee eee eee 


n PITA Ai) Biel Sots deo oe ete 


and 
k ee 
Ud (E 5 = sel 4 2i 2i-1 
ee Pot *n a a 
TP tial 24 2i+1 
mM . 
aft (clave 5 
La one Am Sou omend = a 
eee iet j=l 72541 23 
k s Ley ee ‘ 
a) Seal - aes ee (Uke) 
i=l j=l 723 24-1 


Proof of Theorem 4.3.1 (ii): As noted before, 


g (2) exists iff p LE 
n No vit 


¢3 


Now, for even n, by (4.4.4), 


(De ae eee _ By ys |B ET 
aoe Ge le Mag - 2 a9 1 
V2 182, Vie 12 t i=1 jal eZ 2 j=1 
m—-1 i ans_1 93 
mie See gRL es megane cl 
i=1 = Lee ee 
Since aot] O- Q, therefore 
p_(E(™) > 0 iff u>0oO 
The case when n_ is odd can be dealt with similarly using 
(4.4.7). 0 
Remarks: (i) By noting that 
(a)A5.% (n-1) , . 
p, (Ey) = Pp, (Ey ) \F St a ae S n-L 
(a) \ 2 (n) : 
p,; (E5) = p,; (EB 5) Vd Sai s<en—2 
(4.4.8) 
(i) = (1) poe 


we can calculate all the coordinates of 


Batt Eps Sos Eee and gm) Gabe ad Mh a ae) 


using (4.4.2), (4.4.3), (4.4.4), (4.4.5), (4.4.6) and (4.4.7). 


Cit.) ae ee ba e% OF Eanes’ any exist, must be unique. 


, v oes Aiea 
net i ia ae 
fe a 4 - Ad 
“sae vp, 
ra +9. 
. eat 


: alas 
| ca a 


Sock gph RE ig" 


BY ( iep pae ee 


ie a ae poe 


xe 7 i a 


} j 
f 


| i ene a Oke Opes 


A: Ad bers cept Oa ae 
| Pye a1 re (he. Cop 
Pils Sg i rs ee $e Age ‘ _ - a 5 

. - 7 


| i) > 7 
i oe , . eee ty / 


pi 


Notation: To avoid confusion, we shall denote y (which is defined 


ines 4.3 for (1) ee DY: 2) in the rest of this section and in §4.6. 


Betackss; (i): For 2°< a. <n 


EE exists iff a exists, “ilesse Lif yD) > 0 
(ii) EP a (0, GZ *9 0) 
and 
a 
ee at a Fie ee hak, 
Pai 


always exist. 


Theorem 4.4.1: There exist numbers 


(n) (n) (ny s(n) 
alae eee ilar cg ee 
with 
ig OP Rapin Gh 
such that 
(1) rn <u) see Bee Er ; ges Ee exist 
(2) ne < ym) a Bee 
Lft ae ee eatels Ae exist but not EO? 
) 
(3) ae? < ase? as rae 
ert Ee Eas . hie exist but not ce one 


‘> 
=<? 
a 
ey he 
af Pa a 
ee 
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ee 
- 


“ii ih Wee eve i): “oy 

butideh 22 Hote), arto Lindh ev agian nee been oe Be 
ieee me 

eee RE bod) cers ea ete Sa 7427 oct 5. Fi a Lee +a 


w 2PM. sal 12) 
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a Z yo Suse? 


(i=1) 
; Ga) ea) Ben) , (n) (n) 
Ltt Ey . Ey ; E exist but not Ey ER ueretetars EY 
(n) Stee ym Z yo) 
awk 
: Cr) er) (n) (n) 
tise Eo A Ey exist but not E, ; ; EY 
Proot ; Casevl: n is even, °n = 2m say. Since 
i : : 1 i : ; 
(2m) y ki ea esl = Leis 
One 21 Ome aLoamer ace ae 8 We 6 aa : 
i=l j=l ~2j 2j-1 i=1 RN Sag Lag 
: ‘ k-1 at 
(2k) _ “ 2 82 5n2 224-1 See 
y Pree L 221.0 cae Sf) wut lee ned 
= | 4 J J 2j+1 25 
(for dk m1) 
and 
(2k+1) k 18 k =i 
i, 2j-2 2j-1 dina ’ 
u 7 a6 2 #91 0 ie ER UAne L “24141 0 i eee iui 
=a J ie 2j+1 24 
(Borg. thc) Koc om 1) 
Therefore if we define 
(2m) m i ays > 24-4 m-1 a) Ans 9 24-1 
Eine). 254 9° le Teer ns Pyeaietig lo acrer 
i=k jal ee ae ej! i= j=l = 243 «2j-1 


(ian koe my) 
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: q m—-1 cee: Deas ° 
quiet WD cn eget ky tee ance 
cet | gud hae? 3 925-1 i=k jel Zit 23 
(for Lise <i) 
and 
(2m) 
ei : 
then 
(2m) (2m) (2m) _ 
uy < Uy . ets 2 
and 
on < eee Lt oe) 7 CL ee gee) 


Hence the proof is complete for the case when n is even. 


Case 2: n is odd, n= 2mtl say. “Define 


(2m+1) = = 924-2 (24-1 ee jo 
u mE) 5 hg eras te 85 en a 
2k i=k+1 j=l 724 29-1 i=k+1 j=l 72541 23 
(for 1< k < m-1) 
(2m+1) m i 835 94-1 m Teo ly os 
eee NU . . Soq gel pe Use eae ay ee 
i=k+1 j=l 724 25-1 i=k j=1 725+1 2; 
(for 1 ke m1) 
and 
(2m+1) 
Yom ee 


Then, as is Case 1, we can show that 


(2mt1) (2m+1) (2mt+1) _ 
Wy < Us A eee Pon Q 
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tw 
at 
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oe hi (2m+1) 
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(los te 2m) 
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v4 > i 
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Theorem 4.5.13 


(crt eee 

A 

(n) (n) 

Vv <=> 4 
BRENG 0. Spo, epee 

A C= > Vv <=> A 

Ge5el) (n) (n) (n) 

P, (£3 °°) P,(E, °) p,(E, ) 

(n) (n) (n) 
Py) Py BD eae 

V <=> A <> V 
EECA WAND Eee wean 2 

if n is even, 


and 


Alternating Equilibrium Numbers Phenomenon 


The following two diagrams hold 


ie: 


piod enue Feaetrs 


5 
Q 
Z 
rs 
j 
es 7 
9 
x ; 
~- 
} , 
~® 


pe) = Q 
A 
Py) pM) = 0 
V <=> A 
DAKE ue)) ByGhs aps Cele ye 


A <=> Vv <=> A 


(Gs ae (n) 
BALE, pa ldeg B)(sut meme ske wee) 
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@F542) 
(n) (n) (n) (n) (n) 
Be ae p, (21) Banana: covey, 
A <=> Vv <=> A Ae a V <=> A 
(n) (n) (n) (n) (n) 
Py ae ie) P(E) p,(E ) Peeq ees) PoE ee) 
if im is odd: 
Note: V <=> A or A <=> V indicate that the inequalities are 
equivalent. 
Remark: This phenomenon has been noted by Gard and Hallam [3] for 


n= 4, 


PYOOT ‘6 


f Theorem: It suffices to show 
peo) pen) p,(&.)) Pa Ey) pei”) 
V ZS A <=> V ye Vv << A 
TONS PICO EAC Ce agen ACC) ani a 


if-n is eves, 


' 

é 
a, 
‘a 

ca 


ey 


tr witi 29 anthem: | \ e822 


{ 


athe Lee | 


wer itine ST) saat 


and 
Dey a ka tac ae 
(a5. 4) A <=> v <=> A 
pat "pO Ad (at 


1 


(n) (n) 
Pet nal? ay 
Vv <=> A 


because the remaining equivalences of the two diagrams (4.5.1) and 


morc tollow directly from (4.5.3), (4.5.4) and (4.4.8) 


Case 1: n is even, n= 2m_ say. 


We have to check the equivalence of the following three 


inequalities 
(a) pe”) > ED = 0 
(©) Py, (EQ) > Poy Bony) 
2 (n-1) 
ox En-1 ) 
(1 < k < m1) 
and 
(n) 
(©) Poy By) < Poe Enea? 
(n-1) 
ie ae 


Now by (4.4.4), (a) becomes 


by (4.4.8) 


by (4.4.8) 
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( ear eee: rg i ? ue area 
Siiie? 297-1 21 me, aie co j=l 725 24-1 
m-L a ; 
. 451410 2 2 ~#1 4 > 0 
1 hat pies) 
m ake: 5 F m-1 y OMe: Be i 
ies, ayy > ) ayy g FT : Sit + J Si rieo i 
i=l j=l 725 24-1 i=1 j=l 7244+1 23 


By (4.4.4) and (4.4.5), (b) becomes 


ot 21-2 z Bye grea 
( a . a [a1 oS ) oor 0 Bs a 
Poti =?. 79441 27 i=l fal 724 924-1 
k-1 Zz ays 1 23 
~ Ll aro BG 
i=l j=l 724+1 23 
met ott 0 ae | hae 
> ) we I _£jJ-+ 4). 


i=k 72441 24 j=k+1 725-1 25-2 
By noting that 


1 ea oo Ie = Soon oe 
eS ——— it Go 
1 i 994-1 24-2  i=k "2itl 24 jektl 725-1 25-2 
219. 4) WP Sh 


m-L awe ; 
2j-1 23 

: ADs Op ee , 
i=k qi j=l 225+1 23 


the equivalence of (a) and (b) is then evident. 


On the other hand, by (4.4.2) and (4.4.7), (c) becomes 
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i fr) 
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wy A 
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— io Te a Sa _ 
. +S seer < 
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es eg i 2itt 
f=k “24 21-1. j=k “24 23-1 
k=1 “ale m—-1 isa 
1 © noe 
eet = = i 10 Ps ree a 
See cer eae i=l j=l 725+1 23 
k-1 at ans a=1 
a ) f04 0 : a ] 
i=1 y=], 2424-1 
tom) 2s (Kes 
and by (4.4.2) and (4.4.6), (c) becomes 
m ao. Be 3a ee ae, 
(4.5.6) ) ee SS. 
fed 924 (24-1) j=l “24 24-1 
j paaO Wed 12 Oo) S21 al eo ee 
aia “Ads mle) Otel tae ele 40 
for k= I 
By noting 
m ae f= Pe ee 
1 210 9, _2j 2j+1 
1 ay 994 99-1 iskOAe I Pil gek 2 2 a 
pil 1) 92402141. 
e = 24-2 24-1 
az ) Aoa0 i a 
i=k j=l 725 29-1 


Ine Ges )as and 
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+y a aD C3, 


/ 
‘Nae R tor. 


Ptrd) (on a ‘ Coie. #) bate (SBD § 
| > aa 


80 


m-1 


a al a m-1 eae 
2itl O 24=— Be a e 
a19 ap oes : Dik eer yg il = — 
pele Oe Gee ge oo Gay (el ee Oe 
and 
m Al. . i-l a m ie sy 
2 4 2441 wee 
St | De aoerepeentan ier Be e7 ae erat 
Jee EOE nies) On eat ett $51) 394 24-7 


in (4.5.6), the equivalence of (a) and (c) becomes evident. 


Case 2: Rpts. Odd... n=. 2mt) aesay . 


We ha e to show that the following three inequalities are 


equivalent, 
epee.) Sea eo 
(e) Bae) < Dacre) 
m (n-1) 
=e Oee En ams) by (4.4.8) 
(ls vm) 
Ere ee) Pane, 
Set Bie Sieg at 
(= k sm) 


Since the proof of the equivalence of (d), (e) and (f) is similar to 


that of Case 1 (using (4.4.2)-(4.4.7)), we shall omit the details 


in here. Q 
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84.6 Global Stability 


In this section, we shall give a complete classification of 


the flow of (1): More precisely, we shall show that by decreasing 
(n) 


the value of the parameter H=uU » the equilibria 


5 a 2) 2 6™) 


no i=l | ae ae 


will in turn become globally stable. But before we state and 


prove this theorem, let us consider (I) and (I,) to gain more 
insight into the general case (I). 
Ci) on sans Then (I) becomes 
hs — 
(I,) <a eG LOdea manta, 
Clearly 
Che 
pj CE, #2 =, © 
a 
J: 10 
SHON Ess ee 
ie 


and the phase portrait of (I,) looks like 


(1) (1) 
Eo E 


(1) is globally stable (over Int (R,)) 


Therefore E, 


(4)? n= 2: (I) becomes 


ll 
* 
aN 
ret) 
> 
* 
eh) 


x ) (E £, (x55) 


ae = Xy (ang sr 51%) (= f(x) >*,)) 


he 


1 = 
H 
~ ‘ } a 
‘f i q 
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t ¥ ¥ ij 
ry 4 
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ees ee 
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/ ae 
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~ 5) Se 
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GA aM, a sy Wek 


Clearly, 
2 
E.) =" (0, 0) 
a 
Ee & at , 0) 
17 
and 
2) S20 Wee ak 20> 
ya ion iy pms yas ak 
a a 
Note that ee exists iff _10 » 20 
eg aL 
a a 
Case 1: Bate does not exist, i.e., a ES 
Ae ee | 
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By Coddington and Levinson [35], p. 400, Corollary 1, there 


can be no periodic solution for (I,) iy Int RY 


solution for (I,) in R-. Now the variational matrix of 


+ 


2105-1 e oe Boa eis) 
M = 
51% in he 
C2) 
At Ep 
aso 0 
M = 
0 “An 
(2) 
At E; 
i elo ue 
kh ai 
M= 
A f 2106021 
Zot a 


hence no periodic 
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i’ i ] _ : he 
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, i ; : ; 
Since EC ) is a hyperbolic point whose neighboring flow looks like 
2 
“1 
(2) 
Eo 


then by boundedness of solutions and the Poincare-Bendixson theorem 


(see Coddington and Levinson [34], p. 391), pte? is globally stable 
(over Int (Re). 
a a 
Case 2: EME existsh, 1.e, 5 = > =, Let B(x, »X») ——— 
11 On ep) 
Then 
3(BE,) ; ai 
OX) Xo 
and 
3(BE,) 
ox he 
Z 
so that 
3(Bf,) d(BE,) 
_ lege ee 8G Og. fie (R) 
OX) 9X5 + 


Hence by Dulac's criterion (81.2 (1)) there can be no periodic 


2 
solution of (I,) in wane R.. Moreover, in this case, both oe ) 


and ee are hyperbolic points whose neighboring flows look like 
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X, some eigendirection 
x 
a2) (2) : 
0 uf 
therefore, again by the Poincare-Bendixson theroren, ES is globally 
Stable. 
Theorem 4.6.1: Let a = © and Th =- o, Then 
(n) (n) (st (n) ; 
uy <1 np Re > E+ is globally stable 
COe< k-< a-1) 
Proof: (i) k= n-l. Then gm) exists, let g (9) qh (a Sip he ta) 
a n n 1 n 


and define Vas in $4.2. 


By Lemma 4.2.2, the set 
n 
M = sewey Ln t: Bie ae ees } 


is such that 
ine ak @) outside M 
V.=/0 on. eM 


Therefore to show that x* is globally stable, by the theorem of 


Barbaskin and Krasovskii stated in §1.2 (5), it suffices to show M 
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contains no entire trajectories of (I) except Bean Now suppose 
x(t) = (x, (t), Bl ails x (t)) is a solution of (I) such that 
x(t) eM Werte nO 
Then 
x, (t) = x# 
=> x, (t) = 0 
=> x,(t) = xx by the oe equation so Jl) 
2 Z n 
=> x, (t) = 0 
nd ; 
=> x(t) = xk by thes 2 equation of (I) 
etc. 
(n) 
Hence we must have x(t) = EY Yt. 0 
ae (n) (n) (n) a. 
(ii) Suppose Uy <u <a (OU Kus n=2) <> By 
Theorem 4.4.1, 
(n) (n) (n) : 
Eo - Ey ; ; Et exist 
but not 
(n) (n) 
E 
Eres Seay 
Define the function 
W: Bad naed ss 
where 
n , 
= : k2ac' i <n and 
Bead {x € Rt x, > 0 ( pe ioc nl) 
se, 2 0 (i = i i k+1) } 


ay e oe a est 


seu) dasa (2) Me vbltates See ee eo 
7 + | ‘= a a 5 Cn 


‘S OD ta cal 
He tte 


72 ( £ . Ansa foe ae fog (vw i iw ae 


( : ~ ; ’ Ney 
: i = Oy 


fe ee 


as follows: 


keer 


Wx, S24. = - Saks So 
(xy +++ 5%) a B lx, — Py yyy) - Py Ey) In 
sit 
n 
Piya eee 
i=k+2 7 
Then 
W' (x5. x) 
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Proons KBasy, (c.f.., lemma, 2.2.2 and part (1). 
It then follows from Barbashin-Krasovskii's theorem 
k+l attracts every solution of (I) with initial condition 


in B and hence it is globally stable. QO 


Corallary 4.6.2: Te. <0, then (I) is non-persistent. 


Proof: (By Theorem 4.6.1.)% ~ U 


b)  syvag bine Bl ee Mish 7 
Heroes «! te eve tecd-ereetns SR anti ae ah 
it tOROD ae (g$7af thi 3 7 ; ph Re VOLK iil ipaahen 


Lats obey ah om ae OV 


Fai 
I a0) ; 


nada tagecudcin ts. 2), obi, (eens ee 


yy 


Pa 


aT 


| a 


a: BARS Hi Perot ‘we (i : 


1 


a 
, 
7 
a 
o 
i] 
. 
Ae 
~ ; 
i 
, a 


men it, - 
7 ~ y | 
heidy I a 
ates 
7 a i : a ' = 
- . PR La - a t¢ 2 ? oi : 
j a: eae Up 
a m7 rtd t? . . 


es. 


§4.7 Some Perturbation Theorems 


In this section, we shall consider two perturbed models of 


(I> one with small autonomous perturbation and the other with small 


periodic forcing. 
Model 1: Consider the following perturbed model of (I). 
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where all the specified ai; are >0O, e€ is a small parameter, and, 
f. Geiss. ns fare” Ci. 


We know that if u >0O (where yp is defined in §4.3), 
then (P(0)) has a globally asymptotically stable positive equilibrium 
x*. The following theorem states that under the assumption u > 0, 
then for « #0, an asymptotically stable positive equilibrium point 
Eve) of “(2(e)) bifurcates froms x*. 
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(i) there exists a continuous function «*~> E(e) defined for 
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point of §(P(e)). 
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Gaye Ee) (for + jeih< Eq) is asymptotically stable. 


Proof: Since the proof is standard, we shall only give a sketch of it. 


For (i), define the map 


10 ~ *11*1 7 212% +e £ (xy 5++-5%)) 
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the Jabobian matrix ( —— and @M) issthe variational matrix of 
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(I) ate ox* |. 

(i) then follows from the Implicit Function Theorem (see, 
@.g., Rudin [36], p. 196, Theorem 9.18). For’ (ii) “set sy =x — E(e). 


Clearly (2(Q)) = (I) can be rewritten as 
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x, £, Gx) 
h(x) = 
xf ® 


socgedey oc) =te(yir E(e))i~ s(yet E(O) mats (E(e)),.+ e(B(0)) avsing 
the fact that E(e) is continuous in ¢«, it can be shown that there 


exists a continuous function S(e) of e€ such that 
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On the other hand, by the Mean-Value theorem, it can be shown that 
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Knowing these estimates on h and q, the asymptotic stability of 
E(e) follows immediately from the Variation of Constants formula 


and Gronwall's inequality (c.f., Hale [33], p. 86, Theorem 2.4). JQ 


Model 2: Consider the following perturbed model of (I): 
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where all the specified 75 are >0, ¢€ is a small parameter, 

Pts 1 seen), vare, Cs, sandspertodrcin ete swithe period ‘w. 
The following theorem states that if L270. then for 

e #0, an asymptotically stable periodic solution of (Q(e)) of 


period w bifurcates from the positive equilibrium x* of (Q(0)). 


Theorem 4.7.2: If u> 0, then for small le}, (QCe))» has an 
asymptotically stable solution x(t,e), periodic in t of period w, 
continuous in (t,é), and with x(t,0) = x*. Moreover, there is 


only one such solution for each e. 


BLOOL $ The proof follows trivially from Coddington and Levinson [35], 
p. 348, Theorem 1.1 and p. 350, Theorem 1.2 and the fact that all 
eigenvalues of M have negative real parts (see Corollary 4.3.3 (ii)). 
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CHAPTER V 


EFFECT OF EVOLUTION ON EQUILIBRIUM NUMBERS 


§5.1 Introduction 


Our purpose in this chapter is to investigate the effect of 
selection on the abundance of species of a food chain at equilibrium. 
In particular, we are concerned with the question: "is the driving 
force of selection within a species generally directed so as to make 
the species more abundant, or can selection leave equilibrium numbers 
unaffected or even operate so as to decrease abundance?" 

A method which we now briefly describe for this analysis is 
given in Levins [37]. 

"Suppose that the dynamics of a community of n variables is 
given by the set of n differential equations (1.1.1). There may be 
one or more equilibrium points of the system at which all the f.'s 


are zero. In the neighborhood of any equilibrium point the behavior 


of the system depends on the properties of the community matrix 
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evaluated at the equilibrium point. 
If a genetic variant arises which alters a parameter C. 


it can of course affect many of the species in the system, both directly 


and indirectly. But it will be selected only if 


Therefore the action and direction of Mendelian selection is determined 
by the function f. for the species in question. But the consequences 
of that selection depend on the whole ensemble of species and their 
interrelation, since f. for one species may depend on all the X 5S. 
lf CL is a parameter under genetic control of species i 
(e.g., its bill length, rate of larval development, ability to escape 


predators, etc.), then C, will increase under selection when it enters 


the growth equation 
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for species i with a positive coefficient, that is 210e4 LO 


h 
determine the effect of Cc. on the equilibrium levels of all the 
species in the community, we must differentiate the system of 
expressions on the right of equation (5.1.1) (one such expression for 


each species x.), with respect to Ci. and set the derivatives equal 
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to zero: 
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In §5.2 and §5.3 we shall consider food chains which have an 


equilibrium point at which the community matrix M 
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all the unspecified entries are zero (c.f., Corollary 2.2.2). 
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CL) a= 3: In the following diagram, we have a three-species 


hierarchy of plants x and a carnivore x 


bP herbivore Xo» 3° 


mm 
11 
By Dey 5 
m2 23 


Evolution tewards increased overall viability or fecundity can occur in 
any of the three species. To determine what effeet selection in a species 
will have on its own abundance and that of the other species, we 


evaluate the following determinants. 
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leaves the abundance of the herbivore unchanged; the herbivore has 
more food but suffers heavier predation and merely has a higher 
turnover rate, increasing the abundance of the carnivore; 

(ii) evolution in the herbivore (e.g., towards increased efficiency 
of digestion) does not increase its own population level, and does not 
affect the plant population, but instead increases the abundance of the 
carnivore; 

(iii) evolution in the carnivore affects all three levels: the 


plants increase, the herbivores decrease and the carnivores increase. 


2), a=: In the following diagram, we have a four-species 
hierachy of plants X1> herbivore Xo» carnivore X, and secondary 


Carnivore x 
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n = 3 has been considered by Levins [37]. 
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§5.3 General n 


Lemma 5.3.1: 
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Mtl nt2 


mat? atl 0 


Case 1: n is odd. Then induction hypothesis gives 
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so that’ ‘det A+? = 0. 


Case 2: n is even, n= 2k, say. Then, again by the induction 
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oa, 
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Case, 1: n = 2k+1 is odd. Then 


x k 
det M = m4 (-1) 


by Lemma 5.3.1. 
Gage.2: n = 2k is even. Then 
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where 


1 < Q tf.) gise odd 
(i) when n is odd, det M’* 
J = 0 if j is even 
1 = 0 it: J. tis odd 
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Proof: It follows immediately from Lemma 5.3.2 and Lemma 5.3.3. U 
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abundance of all levels in a manner similar to the cases n= 3,4. In 


fact it is not difficult to show that the follwing tables hold: 
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indicates that selection at the x, species will increase, decrease or 


does not change the abundance of the a species, respectively. 


CHAPTER VI 


EXAMPLES 


One aspect of community ecology that has attracted much 
attention is the relation between the complexity of a food web and its 
Stability. Here increased 'complexity' means more species, more inter- 
action and hence more parameters characterising the interaction; and, 
"stability' means relatively lower levels of population fluctuation, or 
the ability to recover from perturbation, or simply with persistence 
of the system. Evidence that such a relationship exists abound. For 
example, in the study of a marine rocky intertidal community (mussels, 
barnacles and starfish), Paine [23] showed that predators, Pisaster 
ochraceus and Thais emarginata, played the major role in maintaining 
the diversity of prey species. By removing Pisaster, the number of 
species in the marine rocky intertidal community decreased from 15 to 
8 species. However, the issue of whether complexity will promote 
stability is unclear. For example, Elton [38] has drawn together a 
collection of empirical and theoretical observations in support of the 
conclusion that complexity implies srauelity: whereas, Watt [39] gives 
examples of simple natural systems that are stable and of complex ones 
that are not. For a more complete account of the situation and other 
references, see May [15]. 

Our purpose in this chapter is to illustrate the model 


possibility of this question by using Lotka-Volterra ecosystem models 


with 3 and 4 species. 
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§6.1 Example 1 


In this example, we compare the stability of the following 


two ecocommunities: 


2 Cs) 


and 
(A) 
In (A), both predators X» and X, share the same prey x: (B) ds 


the same as (A) except that there is a third predator xy feeding on 
both Xo and x, Suppose (A) and (B) are governed by Lotka-Volterra 


dynamics, then (A) and (B) can be respectively described by the systems 


Dredifterentialvequations (6.1. 1)) and (6.1.2)- 
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where all the specified aij are’ 210% 


Clearly (6.1.1) will have a positive equilibrium only if 


It is known that (see,for example, McGehee and Armstrong [17]) extinction 
occurs for (6.1.1), =1,e., for any choice of oe there is no 
asymptotically stable attractor of (6.1.1) in Int Cae We now show 
that the addition of predator x), has a stabilizing effect, in the 
sense that there is a range of values of the ie Vin factrthigtset 


is a cone containing a non-empty open set of the parameter space 
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such that (6.1.2) has an asymptotically stable positive equilibrium. 


Using the Routh-Hurwicz criterion and with the help of the 


computer, it can be shown that for 


(@.1.3) AAs 


(6.1.2) will have an asymptotically stable positive equilibrium. 
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Tie: fact that the, setof “A’s for which. (6.1.2) has an 
asymptotically stable positive equilibrium is a cone which contains 


a non-empty open set then follows. 
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S052 Example 2 


In this example, we compare the stability of the following 


three ecocommunities: 


(C) (D) (B) 


woere (C)., (D) and (B) are described by (6.2.1). (6.2.2) and (6.172) 


respectively: 
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By Theorem 4.6.1, the Routh-Hurwicz criterion and with the 
help of the computer, it can be shown that for the appropriate 
restriction of A given in (6.1.3) (D) is non-persistent whereas 
there is an asymptotically stable positive equilibrium for (B). This 


indicates that addition of the species x and the two links 


to (D) has a stabilizing effect. 


However, this is not conclusive because 
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both (C) and (D) survive and there is an unstable positive equilibrium 
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§6.3 Example 3 


In this example, we compare the stability of the following 


two ecocommunities: 


C2) = (=) 


(E) (F) 


(6.3.2) x. = x, (an, ot ax =) an axe ) 


(6.3.2) xy = Xy (-ang + 451% 7 ay 3% 4) 


3 = X3(-agq + a3)%) + 859%>) 


Our purpose is to see what happens to the stability of the simple 
food chain (E) (regarding xX, as plant, X, as herbivore and X, 
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Using Theorem 4.6.1, the Routh-Hurwicz criterion and with the 
help of the computer, it can be shown that 


Ci) for 


(E) survives and there is an asymptotically stable positive equilibrium 
for, (E)); 
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(E) survives but there is no positive equilibrium for (F), 
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Gry ant or 


(E) is non-persistent but there is an asymptotically stable positive 


equilibrium for (F). 


Hence we see that depending on the values of the parameters 


240 all possibilities can occur. 


CHAPTER VII 


SUMMARY AND DISCUSSION 


In Chapter I, we proposeda framework (see §1.1) which is 
adequate for the investigation of the dynamics of general food webs. 
The advantage of it is that we can change some of the parameters 


C and see what effect they have on the configuration space. 


poem 
For example, this allows us to study the effect of harvesting, enrich- 
ment, migration or evolution on the ecosystem. Mathematical theorems 
(see §1.2 (1) - (6)) that are useful in the latter chapters are also 
given. 

In Chapter II, ecologically reasonable assumptions (see 
ee Heth) S(yii)sands(1)e= (8)jetare given so thatisygtem, (131.2) 
represents a straight predator-prey food chain with one species 
occupying each trophic level. Even though these assumptions are not 
satisfied by every food chain model studied in literature, they give 
some indication on the adequacy of a model in simulating the dynamics 
of a food chain. Theorem 2.2.1 and its Corollary give conditions under 
which the positive equilibrium of a food chain will be asymptotically 
Stable and for which the variational matrix of (1.1.2) at the equilibrium 
is qualitatively stable. It is to be noted that these results are 
quite general (the conditions only impose restrictions on the signs 
of the entries of the variational matrix and not explicitly on the 
growth functions f\> fos - es tO inu{igli2) jeaerThéicsapplications 


can be found in Chapters III and IV. Theorem 2:orL  isea sufficiency 


theorem on the persistence, non-persistence and survivability of general 
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food webs. Its usefulness depends on the discovery of a persistent, 
extinction or survival function. An application of the theorem to 
Lotka-Volterra food chains, though not given in the above chapters, can 
be found (at least for parts (i) and (ii) of the theorem) in Gard and 
Hallam [3]. There is hope that we can use part (iii) of the theorem 
to improve some persistence results to survivability and that the whole 
theorem is applicable to more general (other than Lotka-Volterra food 
chains) situations, such as food chains with Holling types of predation 
or food webs in general. 

Chapter III deals with intermediate food chains (3.1.1). 
Theorem 3.1.1 states that every solution of an intermediate food chain 
with a carrying capacity in the first species is bounded. A point to 
be noted here is that the method of proof of Theorem 3.1.1 is equally 
applicable to show the boundedness of solutions of a variety of other 
intermediate type food web models (not necessarily food chains). In 
§3.2, the conditions are worked out for the.existence, uniqueness, 
asymptotic stability and instability of the equilibrium points of 
intermediate food chains of length four. Admittedly, the conditions 
obtained (best possible using linearization) are long (hence uninterst- 
ing), nevertheless, the procedure of obtaining them well illustrates the 
use of Routh-Hurwicz criterion, Schwarz matrix theorem and qualitative 


stability matrix criterion (see 51.2 (2), (3), (4)) in model ecosystem 


stability considerations. 


In Chapter IV, we discuss Lotka-Volterra food chains (4.1.1). 
By means of Lemma 4.2.1, system (4.1.1) can be rewritten in a form 
(see (4.2.1)) for which the usual Volterra auxillary function (see the 


Vv function in Lemma 4.2.2) can be constructed. Using this Liapunov 
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function, and under the assumption a > 0, global 


11? 2929 *°*> Aon 


stability for the positive equilibrium of (4.1.1), if it exists, can 
be easily proved (see Theorem 4.2.4). Ina sense, the validity of 


this theorem i : ing' “es 
orem is due to the ‘damping’ terms Ai1> 309° arene 


However, we show in §4.6 that we do not require that much; ari > 0 


is sufficient to ensure the global stability of the positive equilibrium 
of (4.1.1). In 84.3 - 84.6, we concentrate on Lotka-Volterra food 
chains such that Oy ee ee 0 (see (4.3.1)). The results can 


be summarized in the following table: 
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For (ii), we also show that each hypersurface W) (ep O07) (see.844 3) 
contains a recurrent motion. In the case n= 2, it is well-known 
(see Freedman and Waltman [40], for example) that each ve isin fact 
a periodic orbit. Therefore, the following questions arise naturally: 
(a) does Ne contain any periodic orbit in general? 
and 
(b) if so, will each solution of (ii) be periodic (except the 
positive equilibrium itself, of course)? 
For (iii), we know that in general there is no positive equilibriun, 


and that positive equilibria exist when and only when =02) -in «that 


EO 
Case. there will be a half line of equilibria of (iil) in Int (Ry) 


initiating from a point in 
{x € Ry HoesO?, CL = Pewee) Max meen GO) 


Some obvious questions arise: 
(a) can we improve persistence, in case Uy > Oy £0 Survivability? 
and 


(b) what happens to the case Uy = 0; will the system persist? 


For (a), there is hope that Theorem 2.3.1 (iii) is applicable. For (b), 
Freedman and Waltman [6] show that when n = 3 and Wo = 0, then 
every solution of (iii) with initial condition © Lor (Ry) is periodic. 
Therefore we can also ask whether this is true for general n. 

So much for (ii) and (iii). For (i), we have a much 
complete picture (see §4.4 and 84.6). For geometric simplicity, even 
though it holds for general n, we shall only describe our result for 


the case n= 3 here. Consider the Lotka-Volterra food chain with 
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This then answers a question raised by Freedman and Waltman [6]. 

In 84.5, we show that introduction of a new top predator 
into an ecosystem would result in an order interchange of equilibrium 
populations on an alternating scheme. For example, in an aquatic 
System with trophic levels composed of phytoplankton, herbivores and 
primary carnivores, the model conclusion predicts that the addition of 
a secondary carnivore reduces the equilibrium levels of the phytoplank- 
ton and primary carnivore while increasing the herbivore population. 
The two perturbed models (@(e)) and (Q(e)) discussed in §4.7 belong to 
the non-critical case of perturbation theory (since up > 0 implies all 
eigenvalues of the variational matrix of the unperturbed model (I) 
evaluated at the positive equilibrium have negative real parts). There- 
fore standard theorems apply. A (mathematically) more interesting 
problem would be to investigate the perturbed models (P(e)) and 
(Q(e)) when a7 0. Another interesting question is whether the 
perturbed equilibrium E(e) in (@(«)) (we know that it is asymptotic- 
ally stable) is globally stable or not. 

Chapter V is an attempt to investigate what effect does 
selection at each species has on its own abundance and on that of the 
other species in a food chain. Details for food chains of length 
three ona four are given in §5.2. The results indicate that 

(i) Mendelian selection of a species in a food chain may or may 
not increase the abundance of that species, and, 
(ii) the effect of selection on the equilibrium numbers depends 


on the parity (even or odd) of the length of the food chain. 


In Chapter VI, we give some examples to illustrate the model 
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possibility of the widely discussed topic ‘stability verses complexity'. 
The numerical results that we have obtained tend to indicate that 
stability of an ecosystem does not simply depend on the number of 
species present and the number of links between them. ‘Two other 
important factors are: 

(i) the way the species are linked 
and 


(ii) the relative strength of the links. 


However, these examples should not be taken too seriously because 

(i) we have assumed that the addition of extra species or extra 
links does not affect the strength of the original links 
and, 


(ii) our choice of the values for the parameters a has no 


connection with reality. 
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